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Neutrinos may mix with ultralight fermions, which gives flavor oscillations, and with heavier 
fermions, which yields short distance flavor change. I consider the case where both effects are 
present. I show that in the limit where a single oscillation length is experimentally accessible, the 
effects of heavier fermions on neutrino oscillations can generically be accounted for by a simple 
formula containing four parameters, including observable CP violation. I consider the anomalous 
LSND and MiniBooNE results, and show that these can be fit in a model with CP violation and 
two additional sterile neutrinos, one in the mass range between O.f and 20 eV, and the other with 
mass between 33 eV and 40 GeV. I also show that this model can avoid conflict with constraints 
from existing null short baseline experimental results. 
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I. INTRODUCTION 

Since the discovery of neutrino flavor change in a variety of long baseline experiments 1 13 , a new standard picture 
has emerged [TH - f2"4"] . In this picture the neutrino flavor eigenstates e,/i,r are related to the mass eigenstates ^1,2,3 
via a 3-by-3 unitary matrix [25] : 

C12 S12 

-S12 C12 ) I v 2 ) , (1) 
1 

where = cos 6*^, Sij — sin 6^ . Long baseline measurements are consistent with the following values of the three 
angles 6 \f 

tan 2 6>i2 = 0.45 ± 0.05 

sin 2 26 la = 0l° 05 (2) 
sin 2 2023 = LOlg.! 

and neutrino mass squared differences: 



Am 2 2 = (8.0 ± 0.3) • 10~ 5 eV 2 
I Am 2 3 | = (2.5 ± 0.2) ■ 10~ 3 eV 2 

These mass eigenstates are so light and nearly degenerate that in all neutrino experiments neutrinos propagate at 
essentially the speed of light, and the components of the neutrino wavepackets with different mass do not separate 
spatially. The phases of the different mass components oscillate quantum mechanically with different frequencies, 
and, as the different flavors are different superpositions of mass eigenstates, the flavor composition of a neutrino beam 
in vacuum will exhibit spatial variation as a function of L/E, where E is the neutrino energy and L is the distance 
from the source |26j . In propagation through matter the phases are also altered by forward scattering from the weak 
interactions, which alters the flavor change probability in a way which depends on both E and L/E, and differs for 
neutrinos and anti neutrinos [23 [55] . These matter effects are very small in experiments with baseline much shorter 
than 1000 km, unless there exist exotic forces [2"5] - l3"5] . Oscillations are now significantly favored over alternatives such 
as neutrino decoherence or decay [SH EZ] • 

The measurement of the small mixing angle #13 is a primary goal of the current generation of long baseline exper- 
iments. Also sought is evidence for a nonzero value of the CP violating parameter 5, and knowledge of whether the 
pair of states with the smaller mass squared splitting Am 2 2 are heavier or lighter than the third state. 

There have been reports of neutrino flavor change in the short baseline LSND [38 and MiniBooNE [3^MT] experi- 
ments, which would upset this standard picture, as the values of L/E in both these experiments is of order 1 MeV/m, 
which is too small for the small mass squared splittings Amf 2 and A77123 to effect flavor change at the observed 
level. The LSND and MiniBooNE results favor anti electron neutrino appearance in a muon antineutrino beam, at 
different energies and distance but similar values of L/E. The MiniBooNE data on neutrinos disfavors electron neu- 
trino appearance in a muon neutrino beam at the values of L/E explored by LSND, but favors an excess of electron 
neutrinos at higher values of L/E [JO]. With at least 2 additional sterile neutrinos, CP violation in oscillations can 



reconcile the LSND and MiniBooNE anti neutrino results with MiniBooNE neutrino results [4"2T|44| . However recon- 
ciling MiniBooNE and LSND with constraints on muon or electron neutrino disappearance from a variety of other 
short baseline experiments f4*5~W4"5] is more difficult. Attempts to do so have introduced additional exotic ingredients 
beyond neutrino mixing [33 H21 I50H58] . 

Existing studies of neutrino oscillations generally are not sensitive to mass squared differences larger than 1000 eV 2 , 
as the resulting oscillation length is too short to measure. Furthermore mixing with such heavy neutrinos has not 
previously been considered as a mechanism to reconcile LSND and MiniBooNE with short baseline disappearance 
constraints. In this paper we will consider electron neutrino or anti neutrino appearance in a muon neutrino or anti 
neutrino beam in the case where at least one neutrino is so heavy that the associated mass squared difference gives 
rise to an unobservably short oscillation length. 

II. NEUTRINO OSCILLATIONS AND MIXING WITH A NEUTRINO HEAVIER THAN 33 EV 

In this section we will assume neutrino mass differences which are are small enough so that the difference mass 
components of the wavefunction do not separate spatially. In the appendix I review the derivation of neutrino 
oscillation formulae using a wavepacket formalism, and give bounds on the length scale where this assumption to be 
valid. For instance for the MiniBooNE experiment the mass difference must be less than about 20 keV. Assuming a 
unitary mixing matrix, the formula for the probability of electron neutrino appearance in a muon neutrino beam is 
then given by the standard result 



e 2iXil 



i>l 



(4) 



where 

= (mf-m 2 ) L/E 
X v- L - 27 eV 2 m/MeV ' (&j 

Specializing to the short baseline case, where the mass differences among the three light states can be neglected, and 
assuming two additional states, this formula becomes 

Pu^u. = \u e ,u; 4 e- 2 ^ + u e5 u; 5 e- 2 ^ - u e4 u; 4 - u e5 u; 5 \ 2 . (6) 

For anti neutrinos the matrix elements are complex conjugated. This probability can be written as [59j 

Pv^v, = 4|t/ e4 | 2 |^ M 4| 2 sin 2 x 4 i + 4|[/ e5 | 2 |[/ M5 | 2 sin 2 a; 5 i 4- 8 1 C7" eg 1 1 J7 e 4 1 1 ^4 1 1 ^5 1 sinacgi sinx 4 i cos(x 5 i - ar 4i - (j>) (7) 
where 

is a physically observable CP violating phase and for antineutrinos we must replace <p > ~<j). In the limit where the 

5th neutrino is heavy £ 51 varies very rapidly and should be averaged over. The appearance probability then becomes 

Pu^v. MUe4\ 2 \U^ 4 \ 2 sin 2 £4i + 2 1 U e5 \ 2 \ U ^ 1 2 + 4| U e5 1 1 U e4 1 1 U^ 4 \ \ \ sin x ix sin^i + 4>) • (9) 

We can simplify this expression by defining the CP odd quantity /3 

B - 1 tan- 1 ( ™*m*||CUl \ (w) 



and the mixing ratio r 



so that 



and get oscillation probability 



_ \u e5 u; 5 + u ei u; 4 \ 

\u ei u; 4 \ ' [ > 



u e4 u* 4 



|C/ e4 | 2 |C/ AI 4| 2 [2(l-r) 2 + 4rsin 2 /3 + 4rsin 2 (a;4i+/3)] . (13) 



For anti neutrinos we replace (3 — > —j3. Note that CP violation remains observable in the limit of heavy 7715. Note 
also that while with 4 neutrinos the amplitude of oscillations associated with 7714 is 4|[/ e 4| 2 |[/ M 4| 2 , mixing with a heavy 
5th neutrino alters the oscillation amplitude by a factor of r, which may be either larger or smaller. 

In contrast with appearance experiments, in the absence of large matter effects, disappearance experiments are 
much less sensitive to mixing with heavy neutral fermions. Averaging over the short oscillation length associated 
with Amfj, and neglecting the mass differences among the three light neutrinos, the probability for vacuum electron 
neutrino or electron anti neutrino disappearance is 

4|[/ e4 | 2 (l - \U ei \ 2 - |C/ e5 | 2 )sin 2 i 41 + 2\U e5 \ 2 (l ~ \U e5 \ 2 ) (14) 

and the probability for muon neutrino or muon antineutrino disappearance is obtained by replacing e — > fi in the 
preceding formula. Typically, stringent bounds on disappearance are obtained by canceling systematic errors using the 
L/E dependence. Thus U e 5 and are only weakly constrained in the limit where the oscillation length associated 
with Am 2 ,! is too short to give any measurable dependence on L/E. 

The tension between electron antineutrino appearance at LSND and MiniBooNE and short baseline electron and 
muon neutrino disappearance experiments may be reduced by allowing both nonzero (3 and r greater than 1. The ratio 
r is constrained by very short baseline electron neutrino appearance searches. The strongest such constraints come 
from the NOMAD and E776 experiments [BUI EP ■ The NOMAD constraint on u e appearance at L/E < 0.025m/MeV 
implies 

\Ue4,\ 2 \Uf,4\ 2 ((l -rf + 4rsin 2 /3) < 0.0007 . (15) 

For instance for j3 — 0, and |E/ e 4| 2 |E^4| 2 = 2 x 10~ 4 , the maximum allowed value of r is 3. For a 3+1 model with 
S = and r = 1, the maximum probability of electron neutrino appearance as a function of L/E is 8 x 10 -4 . In 
contrast, for r = 3, the probability of electron neutrino appearance maximizes at a much larger 4 x 10~ 3 . Therefore 
the existence of the heavy fifth neutrino makes it is possible to obtain much larger electron neutrino or antineutrino 
appearance probabilities at non zero L/E than would be possible in a 3+1 model. In the next section we will find 
even larger values of r are possible when neutrinos mix with fermions which are too heavy to be produced. 



III. A HEAVY NEUTRINO 



In this section we consider in detail a model with 2 additional sterile neutrinos, one of which is heavy enough that it 
is possible to kinematically distinguish it from the others. In principle it is produced with a reduced phase space, or, if 
sufficiently heavy, it is not produced at all. Mixing with a heavy neutrino will be constrained from loop contributions 
to charged lepton flavor violation, such as /1 — > e~f [55] , but as long as it is lighter than about 40 GeV such a neutrino 
could contribute substantially to the LSND/MiniBooNE anomalies while the charged lepton flavor violation rate will 
be sufficiently suppressed by the Glashow Iliopoulos Maiani [53] mechanism. In the case of neutrino mixing with such 
a heavy neutrino, the beam will not initially be in a pure flavor eigenstate, and the situation can be described in terms 
of a non-unitary mixing matrix for the light states |62. Without the unitarity constraint, CP violation is possible even 
in two neutrino oscillations |64j . In this section we will consider the effect of an additional state which is significantly 
heavier than the other neutrinos, but light enough so that charged lepton constraints on unitarity violation are not 
constraining. 

Neglecting the mass differences among the 3 light eigenstates, but assuming there is a 5th heavy neutrino and a 4th 
neutrino whose mass squared difference with the others is not negligible, the formula for the probability of electron 
neutrino appearance in a muon neutrino beam is 
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giving electron neutrino appearance probability 

Pu^u. = \Ueif\U^\ 2 {(1 -rf +o[(l -r) 2 + 4rsin 2 /3] + 4rsin 2 (x 4 i + /?)} , 



(16) 
(17) 

(18) 



with j3 — > —f3 for antineutrinos. Here a is a phase space factor associated with production of the heavy state, which is 
less than 1, and which is if the state is heavier than the available energy. In this case of a < 1 the constant term is 
smaller than in the previous section, because the contribution from the 5th neutrino to the constant term is reduced. 
Assuming that a = 0, the constraint on r from very short baseline electron neutrino appearance is correspondingly 
weakened, to 



\U ei \ 2 \U^ 



(1 



2r sin 2 /3 



< 0.0007 . 



(19) 



For instance for j3 = 0, and |J7 e 4| 2 |[/ AI 4| 2 = 2 x 1CP 4 , the maximum allowed value of r is 3.8, and the probability of 
electron neutrino appearance maximizes at 8 x 10~ 3 . In contrast, in a 3+1 model with |C/ e 4| 2 |?7 At 4| 2 = 2 x 10~ 4 , and 
no 5th heavy neutrino, the maximum short baseline electron appearance probability would be only 8 x 10~ 4 . 

Also somewhat constraining will be the muon decay rate, which could be affected from and v e mixing with a 
sufficiently heavy neutrino, and which is constrained from lepton universality and from precision electroweak tests[65 . 
For a 5th neutrino which is heavier than the muon, we must require \V e ^\ and \LJ ^ to be smaller than ~ 0.05. 



IV. GENERAL FORMULA FOR ANALYZING NEUTRINO APPEARANCE OSCILLATION 

EXPERIMENTS IN VACUUM 

The results of the previous two sections are easily generalized. Any single neutrino oscillation experiment, is typically 
sensitive to oscillations in a range of L/E which varies by no more than an other of magnitude or so. This suggests 
a simple generalization of the 2 flavor dominance formula. Assuming a single oscillation length is comparable to the 
range of the experiment, much shorter oscillation lengths may be averaged over, and much longer oscillation lengths 
may be neglected, the probability P a ->b of appearance of flavor a in a beam of flavor b in vacuum may generically be 
written 

P a ^ b = sin 2 (2M sin 2 U.27 -J^-^j^ + p \ + .5 cos 2 (20 ab ) sin 2 a , (20) 

where Am 2 is the relevant mass squared difference, ai) is an effective mixing angle, j3 is a CP violating phase difference 
between the two different components (which is allowed to be nonzero if we do not have two flavor unitarity) , and a 
gives the constant term resulting from averaging over short wavelength oscillations, from any nonunitarity, and from 
any neutrinos which are too heavy to participate in oscillations. For antineutrinos the oscillation probability in this 
limit would be 

P- a ^- b = sin 2 (20 afc ) sin 2 (^1.27 ^^^ - p) + .5 cos 2 (20 ab ) sin 2 a . (21) 

This parameterization is chosen to satisfy the generic constraints for CPT conserving vacuum oscillations < P a ->b < 1 
and < (P a _vb) < 1/2. For example, in the five neutrino model of the previous section, we would have 

sin 2 (20 eAl ) = 4\U e5 U; 5 + U eA U; 4 \\U e4 U^\ (22) 
.5cos 2 (26> eM )sin 2 a = |C/ e4 | 2 |C/ M4 | 2 { (1 - rf + a[(l - r) 2 + 4r sin 2 /?] } . (23) 

The CP odd parameter j3 would have the same definition as in that model. Note that a may be as small as in the 
case where r = 1 and a = 0, for arbitrary value of /3. 



V. LSND AND MINIBOONE 



The Liquid Scintillator Neutrino Detector (LSND) experiment at Los Alamos [35] has reported statistically sig- 
nificant (3.8 a) evidence for electron anti neutrinos in a beam produced by the decay of \x + at rest, consistent with 
oscillations of anti muon neutrinos. The MiniBooNE experiment at Fermilab |39H41j which has similar range of L/E 
to LSND, has searched for muon to electron neutrino and anti neutrino appearance. The MiniBooNE electron neutrino 
appearance results showed no excess in the preferred analysis region but do show an excess at lower energies. The 
MiniBooNE anti neutrino data shows an excess which is consistent with a neutrino oscillation interpretation of the 
LSND signal, and which is poorly fit by background. The KARMEN experiment [66] also searched for anti electron 
neutrino appearance in an anti muon neutrino beam, at values of L/E ranging form 0.36 to 0.74, and saw no excess, 
giving a 90% CI limit on the oscillation probability in this region of 0.0017. Several experiments have searched for 
muon to electron neutrino conversion at very short baseline with results consistent with [60] [61] , with the NOMAD 
experiment providing the strongest constraint. 



In order to test whether Eq. (20) could account for the LSND and MiniBooNE results, I have taken the oscillation 
probabilities given in ref. |41] for electron neutrino and antineutrino appearance as a function of L/E to construct a 
X 2 function 

W.„,ffl . V CCTgAmV^) - PTl _ (24) 



where P° xp represents the oscillation probability for bin i extracted from experimental results, P* eory (0 Me , Am 2 , a, (3) 
is given by Eq. (20), averaged over the range of L/E included in bin i, and er^ is the experimental error. I include 
8 bins for LSND, and 9 bins each for MiniBooNE neutrinos, and for MiniBooNE anti neutrinos. I do not use the 



■5 



MiniBooNE data for E < 400 MeV (L/E > 1.37 m/MeV) because of the large systematic error, which should be 
correlated, as I do not have access to the correlation data. Inclusion of these points with the systematic error included 
and treated as uncorrelated makes little difference in the fits, but is not justifiable. I also include in the fit a bin for 
KARMEN, and a bin for NOMAD, with the experimental errors chosen to correspond to the 90% upper bound on 
the average oscillation probability. The total number of fit points included is 28, and there are 4 free parameters. 

The best fit point has Am 2 = 0.40eV 2 , sin 2 (26i Me ) = 0.0083, /3 = -0.123, a = 0, and a total X 2 of 24.14 for 24 
degrees of freedom. A nearly equally good fit may be obtained for any A m 2 in the range from 0.03 eV 2 to 0.60 eV 2 . 
The fit prefers a nonzero value for the CP violating parameter [3. In Fig. |??| I show the lowest x 2 obtainable for a 
given value of Am 2 , with and without the j3 = constraint. The value of \ 2 at 9^ e = a = (no flavor change) is 
45.6. Note that Eq. ( 20 ) appears to give a good fit to all the data on electron neutrino or anti neutrino appearance 



at MiniBooNE and LSND, while being compatible with KARMEN and NOMAD, for a wide range of masses. 

In Fig. [??] I show the LSND and MiniBooNE electron antineutrino appearance probabilities and the MiniBooNE 
electron neutrino appearance probabilities as a function of L/E, together with the curves from 4 points with x 2 < 31. 
Also shown are the constraints from NOMAD and KARMEN. In Fig. |3j I show the region in the Am 2 and 9^ e 
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FIG. 1: The minimum value of the x 2 function described in the text as a function of the Am^ mass squared 
difference in eV 2 . The solid blue line shows the minimum with # Me , a and /3 chosen to minimize x ''■ The dotted 
purple line shows the minimum with # Me , and a chosen to minimize x 2 an d the CP violating parameter j3 set to 0, 
showing that the best fit region has a mass squared difference between 0.03 and 0.60 eV 2 and nonvanishing CP 

violation. 



plane where the x 2 goodness of fit test is within a factor of 10 of the best value ( x 2 l ess than 37.7) for four different 
assumptions about the a and (3 parameters. Because the data has been extracted from the published plots without 
including information about correlations these results should be taken as indicative of the preferred values rather than 
as a definitive constraint region. 

Note that the inclusion of the a and f3 parameters has little effect on the best fit values of Am 2 , but greatly increases 
the preferred region for 9^ e , allowing the effective mixing angle to be much smaller for a given Am 2 than in a 3+1 
model. I do not show the constraints on the allowed region from disappearance only experiments, since, as discussed 
in the previous sections, these depend on U e § and U^, which are only weakly constrained for some values of 7715. I 
conclude that the inclusion of CP violation and nonunitarity in the 3+1 dimensional mixing matrix allows for mixing 
of the 3 active neutrinos with a sub eV mass sterile neutrino to fit all the short baseline electron flavor appearance 
data without necessarily conflicting with muon and electron neutrino disappearance data. Previous attempts |42H44j 
to fit the short baseline data with oscillations among 3 active plus 2 sterile neutrinos were not able to obtain as good 
a fit because only mass squared differences of less than 1000 eV 2 were considered, for which the values of U e $ and 
are more constrained. 
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FIG. 2: Electron flavor appearance probability against L/E in a muon neutrino or antineutrino beam for the 4 
different parameter values indicated, all of which have x 2 < 31 for (24 degrees of freedom). In all cases shown the 

parameter a is zero. The neutrino appearance probabilities are shown in red (med gray) and the anti neutrino 
probabilities in purple (dark gray). The neutrino and antineutrino probabilities differ for the same parameters due 
to CP violation. Also shown are the probabilities extracted from the MiniBooNE neutrino and antineutrino data, 

LSND, KARMEN and NOMAD. 



VI. SUMMARY 



Neutrino oscillation experiments offer an unparalleled window into exotic physics beyond the Standard Model. A 
simple extension of the standard model is to add 'sterile' fermions which are neutral under all gauge interactions. The 
theoretical motivations for such fermions include grand unified theories, Dirac neutrino masses, the seesaw model of 
neutrino mass, supersymmetric models, dark matter theories, and exotic hidden sectors. Such fermions could mix 
with neutrinos, and the mixing angles are not necessarily correlated with neutrino mass. If these exotic fermions are 
light they can appear in neutrino oscillation experiments as a new state, providing a an additional oscillation length. 
Even if they are not light, they can affect neutrino oscillations by allowing the mixing matrix among the light states 
to be nonunitary. 

In this paper I considered the existence of neutral fermions of a wide range of masses, which mix significantly with 
neutrinos, and showed that in the limit of sensitivity to a single oscillation length the usual 2 parameter oscillation 
formula should be generalized to a 4 parameter formula, which accounts for CP violation and for a distance and 
energy independent component to flavor change. This formula has the advantage of covering a wide range of possible 
physics which can affect neutrino flavor change with relatively few parameters. 

I also performed a fit of the electron neutrino and antineutrino short baseline appearance data to the new formula, 
and found a good fit to the anomalous LSND and MiniBooNE results. The best fit parameter region has an additional 
sterile neutrino with a mass squared difference between 0.03 and 0.60 eV 2 , and a nonunitary, CP violating mixing 
matrix, which could result from from mixing with a a second state, with mass between 33 eV and 40 GeV. CP 
violation allows for a good fit to both the MiniBooNE neutrino results and the MiniBooNE and LSND antineutrino 
results. Constructive interference with short distance flavor change term can enhance the amplitude of the oscillatory 
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FIG. 3: The \ 2 < 37.7 region for the Am^ mass squared difference and effective mixing angle is shown in orange 
(medium gray), with a and /3 chosen to minimize \ 2 - Also shown in blue (dark gray) is the preferred region with the 
parameters a and /3 set to 0, which corresponds to a 3+1 neutrino model. The green (light gray) dotted line shows 
the preferred region when /3 = (no CP violation) and a chosen to minimize x 2 i while the yellow (very light gray) 
dashed line shows the preferred region when a = and j3 chosen to minimize \ 2 ■ 



term in appearance experiments, but not in disappearance experiments, allowing for reconciliation of the evidence for 
neutrino flavor change at LSND and MiniBooNE with the lack of evidence from short baseline disappearance searches. 
While this fit is not conclusive evidence for new states or CP violation it is intriguing. 

VII. APPENDIX:STANDARD PICTURE OF NEUTRINO OSCILLATIONS IN VACUUM: WAVE 

PACKET FORMALISM 

Derivations of the standard neutrino oscillation formula have been presented many times. In this section I will use 
a general wavepacket formalism in order to illustrate the approximations necessary for the standard treatment, as 
these approximations break down when mixing with heavier neutral fermions is considered. For simplicity I will only 
consider one spatial dimension. 



Consider a neutrino in a superposition of n mass eigenstates with masses m il where i = 1, 2 . . . n. Its wave function 
may be written 



MO,t)\ 

ip(x,t) = . 

\ip n (x,t)J 

Each component of this wavefunction may be written as a wavepacket, evolving according to 

tl>i(x,t)= [ dp f , 



(25) 



(26) 



where the momentum wavefunctions fi are any square integrable functions of p, Ei (p) — \Jp 2 + m 2 , and we use units 
with h = c = 1. The probability density of detecting flavor b at spacetime point (x,t) is proportional to 



:omm) 



(27) 



where V b * is the amplitude for mass eigenstate i to produce flavor b. In the typical treatment, it is assumed that 
Vbi oc Ubii where U is the unitary matrix which implements the transformation between the mass and flavor eigenstate 
bases. Furthermore, it is assumed that initially 



ip i (x,0) = U ai f(x) 



(28) 



so that initially the neutrino is in a pure flavor a eigenstate. Note that the approximation that the different mass 
components initially have exactly the same wavefunction is not exactly true, since the kinematics will imply that 
different mass components are produced with different momenta distributions. However for the typical case with 
ultralight neutrinos this is a good approximation. 

Although the exact shape of the wavepacket is not important, for simplicity, we take the initial wavepacket to be a 
gaussian which has the particle found in the vicinity of x = at time t = 0. The momentum for each mass component 
centered around pi, with spread 1/d. We take 



f i (p) = N i e- d2 ^ 2 / 2 , 



(29) 



where Ni is a normalization constant, 1/d is the momentum uncertainty and d is the initial position uncertainty. For 
flavor oscillations, d should be small compared with the neutrino oscillation wavelengths 



4p 



(30) 



and the momentum uncertainty 1/d should be negligible compared with the experimental resolution. Kinematics 
of the production will require that pi is a function of i, however any measurement of the momentum with enough 
precision to distinguish two mass eigenstates i and j would require a position uncertainty larger than the corresponding 
neutrino oscillation wavelength A^-. We will assume that for for all i,j 



\Pi-Pj\ < (1/d) < Pi 



(31) 



These bounds are typically extremely well satisfied in realistic experimental situations for the known ultralight neu- 
trinos. 

At a given time t the location Xi(t) of the center of each mass component of the wavepacket moves with group 
velocity Uj, 



dEi 
V{t = ——t 
dp 



p=Pi 



Pi 



d w t 1 - 



The wave packet will also spread out in space as it moves. We may neglect this spread for a time t satisfying 



mft 
2pf 



<C d . 



(32) 



(33) 



In order to measure the momentum with enough precision to distinguish two mass eigenstates would require a position 
uncertainty larger than the ne utr ino oscill ation wavelength. For instance for neutrino energies of order 1 GeV and 
masses of order O.leV/c 2 , Eq. (31) and Eq. (33) give a very weak constraint on the time, t -c 10 22 km/c. On the other 



hand, for oscillations involving 500 MeV neutrinos which mix with a sterile fermion of mass 20 keV, the oscillation 
wavelength would be <~ 0.2^m and, assuming the wavepackets are localized more precisely than this, the different 
mass components of the wavepackets would separate in space after a maximum time of about 500 m/c. Thus even in 
principle oscillations involving a 20 keV sterile neutrino would not occur for 500 MeV neutrinos in the MiniBooNE 
detector. 

For ultralight neutrinos, although the different mass components of the wavepackets may be treated as moving with 
the same group velocity over reasonable timescales, the relative phase between the wavepackets will oscillate with a 
wavelength A, leading to flavor oscillation. 

Specializing to the 2 neutrino case for simplicity of presentation, the probability of finding flavor a at position and 
time (x, t) is proportional to 



P(x,t) 



J dpcos9N 1 e- d2{ P-^ 2 l 2 e l ^ x - E ^ + sm£W 2e -<i 2 (f.-P2) 2 /2 e i( px -E 2 ( P )t) 
For 2 flavor oscillations of ultralight neutrinos we may take N\ = cos BN, N2 = s'mON. Then 

P(x,t) w iV 2 J dpe ip ( x -V (cos 2 Be-^^-P^^e^ 1 

P{x,t) « N 2 J dpe^-^e^e-^tP-P^^t^e 

+ sin 2 6 e - d2 «Pi-P^ 2 P-Pi-Piy 2 )e ~ Z(m ^ m ^ ' ) 

IP1-P2I < (i/<0 



(34) 



(35) 



(36) 



We assume 



(for fixed energy, this is the same as the assumption d <C A), and note that we get very small contribution to the 
integral over p due to destructive interference between different components of the wavefunction unless \p — p\\ < 1/d. 
It is therefore reasonable to approximate 



and get 



P(x,t) w N 2 



dpe 



ip(x-t) 



-d 2 ((p 1 -p 2 )(2p-p 1 -p 2 )/2) 



e-^e-^-fi) 2 ^ ( cos 2, 



sin 



(37) 



(38) 



The condition -C d combined with the fact that the dominant contribution to the integral has \p — p\ \ < 1/d 
allows us to approximate 



e 2p i 



(39) 



This factor then may be taken outside the integral, allowing us to write 



P(x,t) 



P(x,t) 



dpe^^e^e-^tP-P^ 2 / 2 



( cos 2 9 



+ sin 2 Be 2 n ) 



N 1 



dpe 



ip(^-*) P -d 2 (p-pi) 2 /2 



1 - sin 2 (26>) sin 2 



. 2 / (ml - m\)t 



4pi 



(40) 
(41) 



Due to the rapidly oscillating e lp ^ x ^ factor, cancellations between different phase components of this wave function 
will give a small probability of finding a neutrino of cither flavor except when 

\x - t\ < d , 

so oscillations in time will correspond to oscillations in space. Note that have not assumed that that the different 
mass eigenstates have either the same momentum or the same energy expectation value, although we have used a 
small neutrino mass approximation. 
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Recently there has been interest in whether the effects of entanglement of the neutrino with other decay products 
can affect the neutrino oscillation formulae |67H70j . To examine this issue for the 2 flavor case, we consider a muon 
neutrino beam produced from pion decay, which yields a 2 particle state with wave function 

*(*,V,t)=W X > y >V). (42) 
\ip2{x,y,t)J 

Here \tpi(x, y, t)\ 2 is the probability density to find neutrino mass eigenstate i at x and a muon at y. The probability 
to find neutrino flavor a at x and a muon at y is 

| cos 0^i(a;, y,t) + sin 9ip 2 (x,y,t)\ 2 . (43) 

If ip(x,y,t) does not factorize into x{Vi t)^( x t t) we sav the muon and the neutrino are in an entangled state. To find 
the probability of finding neutrino with flavor a at x, independent of where the muon is, we integrate over y to get 

P(x,t) = J dy\cose^ l {x,y,t) + sin6il) 2 {x, yi t)\ 2 . (44) 

Neglecting interactions, each of the wavefunctions ip^ satisfies the Schrodinger equation for 2 free relativistic particles, 
that is a general solution has 

*i{x,V,t) = JdqJ dpfi&qy^y-^+^t (45) 

where 



E(q) = ^q 2 + M 2 , (46) 

and E does not require an i index because the muon only comes in 1 mass eigenstate. The fact that the expectation 
value of the momentum of the muon depends on which neutrino eigenstate is produced is reflected in /1 7^ f 2 . This 
is the most general possible solution. To get a wave packet for each particle, with both localized in the vicinity of 
the origin at t = 0, we choose fi{p,q) to be square integrable functions which do not contain any rapid change of 
phase as p, q are varied. A product of gaussians or any similarly shaped function centered around (jPi,qi) will do, with 
width in both directions greater than 1/A, consistent with position uncertainty less than A. As before, to get neutrino 

oscillations, we assume CrfcA, and \pi — p 2 \ ~ \qi — q 2 \ ~ 0(1/A) ~ 0(\m\ — m 2 \/pi). 
The probability density of finding a muon flavor neutrino at x, t is 



P(x,t) = JJJJJ dydqdq'dpdp' (cos Of i(p,q)e- iEl ^ t + sm0f 2 (p,q)e- lE2{p)t )e l{xp+y<! ~ E( - q)t) 

x(cosef*(p',q')e lEl< - p ' )t +smef;(p\q')e lE2 ^ t e- l{xp ' +yq '- E(q ' )t) ) . (47) 

We can do the y integral, which will give a delta function S(q — q'). We then can do the q' integral. Note that when 
we do this, the dependence on the energy of the muon goes away, and we are left with 

P(x,t) = JJJdqdpdp' (cos 9 f^p, q)e- iEl( - p}t + sm9 f 2 (p,q)e^ E2( - p '> t )(cose f*(p' ,q)e iEl ^'^ (48) 

+ sinef;(p',q)e lE ^ p "^ t )e lx( - p - p " > ) 

- JJJdqdpdp'( cos 2 e(h(j>, q)h{p\ qfe^-^-^^-^'W (49) 

+ sin 2 0/ a (p, q)f 2 (p', q)*e^ p ~ p ' ) ~ (E ^ p) ~ E2{p ' ))t) (50) 
+ sin(20)5R(/ 1 (p,g)/*(p', g )e"^')-( Bl ^- B2 ^')) 4 )) 

W JJJdqdpdp' (COS 2 e(h(p, q)h{p', ? )*e i ((- t )(P-f')+(m?(l/(2p)-l/(2p')))t) 

+ sin 2 6f 2 (p, q)f 2 (p', ? )* e i((»-t)(p-p')+(™l(V(2p)-i/(2p')))t) (51) 
+ sin(26')K(/ 1 (39,g)/*(p', g )e l ( :l: - t )^ p ') + ( m i/( 2p ^" l ^( 2p '» t )) . 

Note that as before the location of the expected position in space is correlated with the time, due to an oscillating 
phase factor giving destructive interference in the integral unless x ps t ± d/2. We can therefore consider flavor 
oscillations in time, and perform an integral over x. Since the center of the wavepacket moves with constant velocity, 
oscillations in time will correspond to observation of oscillations in space. 
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We can find the probability as a function of t to detect a neutrino of flavor a 

/oo 
dxP{x,t) . (52) 
-OO 

The integral over x will provide a delta function 5(p — p') and so we can also do the integral over p' . The result is 

P(t) ^JJdqdp(co S 2 9\h( P ,q)\ 2 + sin 2 9\f 2 (p,q)\ 2 + sin(2^)5R(/ 1 (p, g)/ 2 *(p, g)^^"^)*/^)) . (53) 

Note that the last term, which represents the interference, is now the only term which oscillates in time and it oscillates 
at the canonical rate. Its size, which gives the amplitude of the oscillations, may be reduced by the entanglement, 
which reduces the overlap of the momentum space wavefunctions for different mass components. For sufficiently light 
neutrinos the difference in the central values of the momenta in the momentum space wavefunctions will be much less 
than the intrinsic uncertainty in the parent particle energy due to its finite lifetime and the effect of the entanglement 
can be neglected. 

For sufficiently large mass differences the interference terms between different mass components of the neutrino 
wavefunction will be suppressed due to the near orthogonality of these components, and the amplitude of oscillations 
will be suppressed, without affecting the average amount of flavor change. 
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